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Abstract 

We study the Zeeman effect on entanglement of non-equilibrium finite-spin systems with exter¬ 
nal fields using a method based on thermofield dynamics (TFD). For this purpose, the extended 
density matrices and extended entanglement entropies of two systems with either non-competing 
or competing external fields are calculated according to the dissipative von Neumann equation, 
and the numerical results are compared. Consequently, through the “twin-peaks” oscillations of 
the quantum entanglement, we have illustrated the Zeeman effect on the entanglement of non¬ 
equilibrium finite-spin systems with competing external fields in the TFD algorithm. 


PACS numbers: 03.65.Ud, ll.lO.-z, 05.70.Ln, 05.30.-d 


1 Introduction 

The entanglement of quantum states is a correlation between multiple systems that is peculiar to 
quantum mechanics [T] . The behavior of entangled quantum states is called “quantum entanglement”. 
It plays important roles in quantum computation and quantum information [5], and is useful in 
applications of the AdS/CFT correspondence [3J|3]. Entanglement entropy is a measure of the strength 
of quantum entanglement. It has been used as an order parameter of quantum spin systems [SJ E]. 

A different method to analyze quantum entanglement using thermofield dynamics (TFD) [7l |8l |9] 
was proposed in Ref. [10], and its applicability confirmed in Ref. In this treatment of quantum 
entanglement with TFD, an extended density matrix is defined on the doubled Hilbert space (physical 
and ancillary Hilbert spaces), and examined for some simple cases (TU]. The TFD-based method allows 
the entanglement states to be easily understood, because the intrinsic elements caused by quantum 
entanglement can be extracted from the extended density matrix in this formulation. Consequently, it 
was found that the intrinsic quantum entanglement can be distinguished from the thermal fluctuations 
included in the definition of ordinary quantum entanglement at finite temperatures. Based on the 
analysis presented in Ref. uni, it was argued that the general TFD formulation of the extended density 
matrix is applicable, not only to equilibrium states, but also to non-equilibrium states. The extended 
density matrix was calculated as a simple example in Ref. nni for the case of equilibrium finite-spin 
systems with external fields. In Ref. [m, it was shown that the extended entanglement entropy, which 
was obtained by using the extended density matrix, could be decomposed into parts from thermal 
(but classical) fluctuations and from quantum entanglement, as in Eqs. (fTBl) ^ (fT51) . In the present 
communication, we prove that the value of the extended entanglement entropies of the non-equilibrium 
finite-spin systems is positive semi-definite. We then obtain the extended entanglement entropies of 
TFD as a well-defined measure of quantum entanglement. 

Moreover, in Ref. [1^, the authors gave an example of a frustration effect on the entanglement of 
equilibrium finite-spin systems that is caused by competing external fields. They compared a system 
that does not contain a competing effect between the interaction and the external fields with one that 
does contain it. They concluded that this competing effect is a kind of frustration, in relation to 
the partial recovery of the broken symmetry of the spin inversion, by using the equilibrium density 
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matrix. However, in Ref. |10] . because only a two-spin system is considered, we feel that it is incorrect 
to refer the effect of external fields as “frustration”. In addition, the effect of the external fields was 
analyzed only in an equilibrium two-spin system, and not in a non-equilibrium two-spin system m- 
The effects of external magnetic fields have generally been considered among the Zeeman effects. 
Thus, the Zeeman effect on the entanglement of non-equilibrium systems with external fields remains 
of interest. In the present communication, we therefore investigate, exhaustively, the extended density 
matrices and entanglement entropies of non-equilibrium spin systems with both non-competing and 
competing external fields based on the general TFD algorithm. By comparing the quantities for the 
non-competing and competing cases, we demonstrate the Zeeman effect on the dissipative dynamics of 
the entanglement through “twin-peaks” oscillations of the extended density matrices and the extended 
entanglement entropies. Furthermore, the origin and the generality of the “twin-peaks” oscillations 
are argued by using the extended entanglement entropies. 

This paper is organized as follows: we introduce the extended density matrices of the non¬ 
equilibrium systems with external fields in the next section and examine their properties. In section 
3, we obtain the extended entanglement entropies of non-equilibrium spin systems with competing 
and non-competing external fields, and discuss the numerical results. The last section is devoted to 
discussion and conclusions. 

2 Extended density matrices of non-equilibrium finite-spin 
systems with external fields 

The Zeeman Effects was originally studied in the Ising spin-glass model, which is described by the 
Hamiltonian 


n = - 


(hi) 


■ S, 


hi ■ Si, 


( 1 ) 


where Jij is the strength of spin interaction, hi is the external magnetic field, and (i, j) express the 
pairs of nearest neighbors. Here, the Hamiltonian is the most simple separation of Eq. O- Consider 
the S' = 1/2 spin system described by the Hamiltonian [TU] 

T-L := —JSa ■ Sb — gfJ-s (HaS^ + HbSb) , (2) 

which involves the spin operators, Sa = {5%, S^, S^) and Sb = {S^, S^, S^) of the subsystems A 
and B, respectively, where g is the Lande factor, hb the Bohr magneton, and Ha and Hb are the 
components of the external fields conjugate to Sa and Sb, respectively. To examine the Zeeman 
effect on the entanglement. Ref. m recognized that the two systems with Hamiltonians, and 
should be compared for the cases of a non-competing external field Ha = Hb — H and competing 
external field Ha = —Hb = H, respectively. As can also be seen from Eq. ([2]), the spin inversion 
symmetry of these systems is broken by these external fields in the Hamiltonian. In a forthcoming 
analysis exploring the Zeeman effect on the dissipative dynamics of entanglement, the partial recovery 
of this broken symmetry will play an important role. 

The state, |s), of the total system is defined by the direct product, |s) = |sa, sb) = |sa)|sb)- Using 
the base {| -I- -|-), | H—), | —h), |-)}, the matrix form of the Hamiltonian ([5]) is expressed as 

= -JSa ■ Sb - ggsH {SI + 

= - ggBH^ I + +)(+ + I + + J^bH^ \ -)(-| 

+ ^(l+ -)(+-l + |-+)(- + l)-^(l + -)(- + l + l-+)(+-!) (3) 
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for the non-competing case and 


= -JSa ■ Sb - g^isH {SI - S^) 

= (I + +)(+ + 1 + 1-)(-I) 

+ “ ggsH^ I + -)(+ “ I + + gg-sH^ \ - +)(- + I 

-^(l + -)(- + l + |-+)(+-|) (4) 

for the competing case. 

Next, consider non-equilibrium systems with dissipation, described by the Hamiltonian of Eq. ©■ 
The time dependence of the ordinary density matrix, p“(t), of these systems is given by the dissipative 
von Neumann equation HHS] 

= [H“,p“(t)] - e (p“(t) - , (5) 

where e is a dissipation parameter, a “nc” for the non-competing case or “c” for the competing case, 
and p“q the ordinary density matrix of the equilibrium systems m- The solution of Eq. is then 
expressed as 

p-{t) = e-^‘C/“t(t)poC/“(t) + (1 - (6) 

for an arbitrary initial density matrix, po: where the unitary operator, U°‘{t) := denotes 


[/”(() = + + +)(+ + 
+ exp I--)(-- 

+ coey(| + -)(+-| + |-+)(- + |) 

-iein^|(|-+)(+-| + | + -)(- + |) 

for the non-competing case and 

I 

e—‘/2(|++)(+ + ! +I--)(--!) 


u\t) = 

/ 
V 


-I- 


cos 


V 


(+uj‘^h? 


2h 




x|+-)(+- 

/ 


V 


(ty'4g^H^fil+uj^hA I sm 2 fi j 


x|-+)(- + | 


(|-+)(+-| + l + -)(- + 


(7) 


( 8 ) 
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for the competing case. Here, w := J/h. The explicit expression for /o“(t) in Eq. ([6]) is complicated 
for an arbitrary initial condition and it is difficult to understand its physical meaning. Hence, we will, 
for the sake of simplicity, hereafter confine the discussion to the initial condition po = \ -\ —)(H— |. 
Inserting Eq. © or ®, along with the initial condition, into Eq. ®, we then obtain 


"“W = ? ++><+ + I + I - -X- - I) 

/ (e'^* — 1) cosh (y) 


^cosh (y) + cosh(/i)e'^/2 

^ — 1) cosh (y) 

^cosh (y) + cosh(/i)e'^/2 
' e'* (e"* - 1) (e^ - 1) 

gif ^g2/i _|_ gh _|_ X) -(- 


+ cosuit + 1 I I H—)(H— 


- coswt + 1 ] I —!-)(- + 


— isinwt ) I H—)(—h 


/ - 1 ) (e^ - l) V \ 

1 — u f Oh - h -- h * sinwt 1 - — 1 ) 

y + 1) + J j 


for the non-competing case or 


XW ( e»4cosh(L) (I + +><+ + I + I - -)<- - I) 

16Le^+^*/^ sinh( f )(/i sinh(L)-|-L cosh(L)) 1^2 ( 2U\ _ If 2 , o r 2 

2eK/2+L f g2L f i *^0® I if j ^ -Hon 


4L2 


^ I ”1 )(“l I 

' sinh(f )(Lcosh(L)-?isinh(L)) ^ ^^2 „■ 2 t Lt\ 

2 gK/ 2 +I. , e 2 I,+i I- n. sm ^ j 


2L2 

Kf4:hsm^(^) 2 - 1 ) - 1) . f2Lt 

L + 2e^/2+i-b e2i-h 1 

X I + -)(- + I 

KfAhsm^(^) 2 - 1 ) - 1) , f2Lt 

"^4l( L + 2e^/2+i-b e2i1 


-+)(- + 


(9) 


x|-+)(+-| , 


( 10 ) 


for the competing case, where the scaled time, cut, and dissipation rate, e/ui, have been used. In 
addition we have introduced the definitions K := (3J = (iuih, h := PgpsH, and L := jA. 

The extended density matrix, /3“, in the TED doubled Hilbert space was defined in Ref. [10] as 

:= \^){n 1^) := (11) 

S S 

using the ordinary density matrix, p“(t) in Eq. (|n|). Here, {|s)} is an orthogonal complete set in the 
original Hilbert space and {|s)} the same set in the ancillary Hilbert space of TED [T31[T5]. If the 
entangled subsystems A and B are being examined, each of the |s) and |s) states is represented as a 
direct product: |sa,sb) = |sa)|sb) and |sa,Sb) = |sa)|sb), respectively. 

According to the matrix algebra of p“(t) developed in Ref. (TU], we then obtain the extended 


4 



















density matrix, 


Pl = b2i\+){+\\+){+\+b^2\-){-\\-){-\ 

+ &S(l+)(-ll+)(-l + l-)(+ll-)(+l) 

+ &qe(l+)(+l|-)(-| + |-)(-|l+)(+l), (12) 

where the matrix elements, ^cfi ^'^‘1 ^qei obtained as analytic functions of H, e, t, and /3, 

respectively. They correspond to the two diagonal elements (dl and d2), classical and thermal fluctu¬ 
ations (cf), and quantum entanglement (qe) of respectively. Their expressions are so complicated 
that we show only the numerical results in Fig. [T] for the non-competing cases and in Fig. [5] for the 
competing cases. The asymptotic behaviors of &di,&d 2 >^cfj ^^*1 ^qe at t —>■ oo correspond to those of 
the equilibrium systems as follows: 


lim b 

t—^OC 


nc _ 

dl — 


lim 6^2 = 

t—¥oo 


lim = 

t—¥CO 


4e^+j^ + (e^/2 + l)^ 


4e^(2 cosh(/i) + 1) -|- 4’ 
(e^ + 4) -I- 

4 ((e^ -I- + 1) e-^ + e^) 

(e^ + 1 ) (e^/^ + 1 ) 

2 ((e^ -I- + 1) + e^) 


(e^/ 2 - 1 )" 

and lim b^l = ^ -- 

t-yoo a 4e^(2 cosh(/i) + 1)-I-4 


(13) 


and 


lim ^di = 


t—^OO 


32 L 2 + cosh(L)) 


(4L((e^^ - l) \/4L2 _ ^2 


lim 6 d 2 = 


lim blf = 


+ 2L (2e^+^ -f -1-1)) - (e^ - l)^) 

T 

■(4l(2L 


2e^+^-b 


32 L 2 + cosh(L)) 

- - 1 ) v' 4L2 _ j ^ 

(gi + 1) eK^+2i) 


-b 1 


2 e^ 


+L 


and lim = 

t^oo 


+ + 1 
(e^ - 1 )^ 


16L2 


2ef-b -b 1 


(14) 


respectively. The asymptotic values with K ^ = {fiJ) ^ = 0.7 and h = fUgfisH = 3/7, which are 
obtained from Eqs. (US and (HI, are shown in Figs. 1 and 2. It is worth mentioning that the findings 
of Figs. 1(a), (b), and (c), and Figs. 2(a), (b), and (c) may be in harmony with the relaxation processes 
of entanglement by dissipation. 

As can be seen from Fig. 1(d), which shows the non-dissipation case (e = 0), vanishes identically. 
On the other hand, b^'^, 6 ^( 2 , and show a kind of “classical-quantum crossover” oscillation in 
Ref. [TU]. However, as can be seen from Fig. 2(d), which shows the non-dissipation case (e = 0), 
vanishes identically, while the b^^ curve displays a “twin-peaks” oscillation, which is a new type. 

In Ref. [TU], the following results were obtained, illustrating the equilibrium case: i) In the non¬ 
competing systems (for zero temperature), the external field, H, breaks the spin inversion symmetry 
and the parameter 6 ))“ = 1, because the quantum entanglement 6 “ = 0. ii) In the competing system, 
the level splitting creates a non-zero finite entanglement, 6 “^ 7 b 0 , for zero temperature, even in a finite 
external field, H. hi) The parameter, , which expresses the probability weight of the up-state, 
is smaller than the maximum value, I.O. According to these observations, Ref. m demonstrated a 
typical example of entanglement caused by the external field. 
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Therefore, also in the present analysis it seems reasonable to conclnde that the competing “twin- 
peaks” oscillatory behavior is a consequence of the partial recovery of the spin inversion symmetry 
even in the dissipative dynamics. 



0 2 4 6 8 10 12 14 Cdt 0 2 4 6 8 10 12 14 tji 

(a) (b) 



0 2 4 6 8 10 12 14 Cdt 0 2 4 6 8 10 12 14 tjt 


(c) (d) 

Figure 1: Time dependence of the matrix elements, and fcqe(-•-), in dis¬ 
sipative and non-dissipative systems with K~^ = (/3J)“^ = 0.7 and h = = 3/7, for the 

non-competing cases. Parts (a), (b), and (c) show cases with a scaled dissipation rate e/w = 1, 0.1, 

and 0.01, respectively. The dot-dot-dashed lines (-) in parts (a) and (b) represent the asymptotes of 

the &di, ^qe curves, respectively. In part (d), which is the non-dissipation case (e = 0), 

vanishes identically. On the other hand, s-nd show a kind of “classical-quantum crossover” 

oscillation. 
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0 2 4 6 8 10 12 14 CJt 0 2 4 6 8 10 12 14 Oli 


(c) (d) 

Figure 2: Time dependence of the matrix elements, &dl(—),&d2(-).&cf(-), and 6qg(—), in dissipa¬ 
tive and non-dissipative systems with K~^ = = 0.7 and h = = 3/7 for the competing 

cases. Parts (a), (b), and (c) show cases with a scaled dissipation rate of e/uj = 1, 0.1, and 0.01, 

respectively. The dot-dot-dashed lines (-) in parts (a) and (b) represent the asymptotes of the 

5gf, &J 2 ) ^'^d foqg curves, respectively. In part (d), the non-dissipation case (e = 0), vanishes 
identically, while the 6qg curve displays a “twin-peaks” oscillation. 
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3 Extended entanglement entropies of non-equilibrium finite- 
spin systems with external fields 

The extended entanglement entropy is defined as m 

5“:=-fcBTrA[pllogpX], (15) 

using Pa in Eq. (IT^ . Inserting Eq. (IT^ into Eq. (ITH1) and subsequent simplification eventually yields 


QCt _ OCX I QCX 

~ *^01 ‘^qe’ 


(16) 


where 


Sci ■= -kB I V 4 + ( 6(^1 - arccoth 


4- h 

^dl ^ ^d2 




ua I Lo 
^dl ^ ^d2 


log (&dl&d 2 - (&cf)' 


(17) 


and 


:=-2fcB6“3log6 


cx 

qe5 


(18) 


respectively. In Eqs. (HU), (HZD, and (HU, the expressions for the classical and thermal fluctuation 
parts, and the quantum entanglement part, S'qg, also involve analytic functions of t,/3, e, and 
uj. However, these expressions are quite complicated. Therefore, we show the numerical behavior of 
5“g, and 6 “g for a few scenarios with K~^ = (/3J)“^ = 0.7 and h = = 3/7 in Fig. [3]for 

the non-competing case and in Fig. H] for the competing case. The asymptotic behaviors of and 
5qg as t —>■ oo correspond to those of equilibrium systems obtained by inserting Eqs. o and H3 
into Eqs. (ITOl) . (fT71) . and (fTSl) . The asymptotic values of S°‘ and 5qg with K~^ = (/3J)“^ = 0.7 and 
h = PgfiBH = 3/7, obtained from Eqs. (fTHll . (fT71) . and m, are displayed in Figs. 3 and 4. The data 
in Figs. 3(a), (b), and (c), and Figs. 4(a), (b), and (c) are in good agreement with the relaxation 
phenomena of entanglement by dissipation. 

In the non-dissipative systems, (e = 0), in Eq. (IT^ reduces to 


'^qe = sin^ t ■ log (4 csc^ t) , 


for the non-competing case and 


, fa(coq,/4FTTt) + 8tf + i) ^._^, / 

Uh^ + lf V2 J 


X log 


(4/i2 -P 1)" 

2csc^ (^\/4/i2 -I- It) 
cos (\/4/i2 -p it^ -p 8 /i 2 _|_ I 


-h 2 log(4/i2-p 1) 


(19) 


( 20 ) 


for the competing case. As can be seen from Eqs. (ITOl) and (001) . both the K and H dependences of 
Sqg at e = 0 disappear and the K dependence of 5qg at e = 0 disappears. The time dependence of 
and ,5qg at e = 0 is shown in Fig. [SKd) for the non-competing case and in Fig. |4Dd) for the competing 
case. It is apparent from these figures that the curves (5“ and S^^), showing entanglement, have the 
same phase. However, their amplitudes differ. Especially, the H dependences of 5qg at e = 0 are 
displayed in Fig. 5. From this observation we infer that the “twin-peaks” oscillation is an example of 
quantum entanglement caused by level splitting. 

Next, we consider the concurrence of non-equilibrium finite-spin systems with external fields. The 
ordinary entanglement entropy. S'", is defined by 


S“ := 


feuTrA [pa log Pa] > and p% := TrBp“(t). 


( 21 ) 












The insertion of Eq. m into Eq. (EH yields 



where C°‘ is the concurrence |16] for the competing or non-competing cases. The time dependence of 
C°‘ is displayed in Fig. |2] for the non-competitive case and in Fig. H] for the competitive case (in units 
of fee = 1). As can clearly be seen from Figs. |3] and 01 C“ includes not only the contribution from the 
quantum entanglement, but also the contributions from the classical and thermal fluctuations in the 
non-equilibrium case. However, this fact is not manifest in the above expressions for 5"“ and (7“. 






Figure 3: Time dependence of -), «Sqe(—), and C“(*") in dissipative and non-dissipative 

systems with K~^ = {j3J)~^ = 0.7 and h = = 3/7 for the non-competing cases. Parts (a), 

(b), and (c) show cases with scaled dissipation rates of e/w = 1, 0.1, and 0.01, respectively. The 

dot-dot-dashed lines (-) in parts (a), (b), and (c) represent the asymptotes of the and 

curves. In part (d), which is the non-dissipation case (e = 0), all the curves have the same phase. In 
these and all subsequent plots, A:b = 1 - 
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Figure 4: Time dependence of S^{ -), 5qg(—), and in dissipative and non-dissipative systems 

with K~^ = = 0.7 and h = = 3/7 for the competing cases. Parts (a), (b), and (c) 

show cases with scaled dissipation rates of e/w = 1, 0.1, and 0.01, respectively. The dot-dot-dashed 

lines (-) in parts (a), (b), and (c) represent the asymptotes of the S'^ and curves, respectively. 

In part (d), which is the non-dissipation case (e = 0), all the curves have the same phase and show a 
kind of “twin-peaks” oscillation. 
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4 Discussion and Conclusions 


In this communication, we have examined the extended density matrix and the extended entanglement 
entropies of non-equilibrium spin systems with non-competing and competing external fields based 
on the TFD formulation. These results are summarized in Figs. [I] ^ [S] As is evident from the 
results, the values of are positive semi-definite. In particular, according to Eqs. (12) and (13), 
the equilibrium-extended entropies, lim 5'“, vanish at ^ 0. From these observation we infer 

t —)-00 ^ 

that S'qg are well-defined measures of the quantum entanglement. This is also the condition that the 
entanglement measures of general quantum systems should satisfy. 

Figures 1 and 3 show that the behaviors of 5“, and 6“ differ little from those with no external 
fields [101 ITT| . and are not indicative of the Zeeman effect. This is due to the lack of competition 
between the interaction, —JSa ■ Sb, and the external fields, —g^sH {S^ + S^), in the non-competing 
Hamiltonian, 'H™ in Eq. ®. The spin inversion symmetry is still unbroken for the non-competing 
case. On the other hand, it is clear from Figs. 2 and 4 that the behaviors of S^, and differ 
from those with no external field. Specifically, Figs. 2(d) and 4(d) show a “twin-peaks” oscillation 
with the same phase. This “twin-peaks” oscillation may be caused by the H dependences of at 
e = 0, which are shown in Fig. [S] These results show that the competing Hamiltonian, in Eq. ®, 
contains competing effects in a non-equilibrium system leading to the partial recovery of the spin 
inversion symmetry in the dissipative dynamics. This is recognized as the origin of the “twin-peaks” 
oscillation in the present model. 

In the present communication, the Hamiltonian was the most simple separation of Eq. O- Using 
the replica trick and the Sherrington and Kirkpatrick mean-field solution nzmiiin], the free-energy 
and the order parameter of the equilibrium system of Eq. © have been examined. However, these 
methods involve several serious problems, such as replica-symmetry breaking, ergodicity breaking, 
negative entropy at low temperatures, etc. On the other hand, the method presented here should 
be independent of the problem. Therefore, one can extend the method developed here to study the 
quantum entanglement in other spin systems. Furthermore, to clarify how general the appearance of 
“twin-peaks ” oscillations is in spin systems, it is necessary to try the present method based on TFD 
in the above general spin-glass models. These extensions are currently under consideration. 
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